About coherent structures in random shell models for passive scalar advection 
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A study of anomalous scaling in models of passive scalar advection in terms of singular coherent 
structures is proposed. The stochastic dynamical system considered is a shell model reformulation 
of Kraichnan model. We extend the method introduced in JlJ to the calculation of self-similar 
instantons and we show how such objects, being the most singular events, are appropriate to capture 
asymptotic scaling properties of the scalar field. Preliminary results concerning the statistical weight 
of fluctuations around these optimal configurations are also presented. 
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Experimental and numerical investigation of passive 
scalars advected by turbulent flows have shown that pas- 
sive scalar structure functions, T2 n (r) have an anomalous 
power law behaviour: T 2 „(r) = ({8{x + r) — 9{x)) 2n ) — 
{(S r 9(x)) 2n ) ~ r C(2™) 5 wnere f or anomalous scaling we 
mean that the exponents C(2n) do not follow the di- 
mensional estimate £(2n) = n((2). A great theoretical 
challenge is to develop a theory which allows a system- 
atic calculation of £(n) from the Navier-Stokes equations. 
Recently |] , it has been realized that intermittent power 
laws are also present in a model of passive scalar advected 
by stochastic velocity fields, for n > 1 The model, 

introduced by Kraichnan, is defined by the standard ad- 
vection equation: 
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u • 89 = nA9 



(1) 



where u is a Gaussian, isotropic, white-in-time stochas- 
tic d-dimcnsional field with a scaling second order struc- 
ture function: ((ui(x) — Ui(x + r))(uj(x) — uj(x + r))) = 
Dor €((d + £ — l)Sij — ^Tirj/r 2 ). The physical range for 
the scaling parameter of the velocity field is < £ < 2, <f> 
is an external forcing and n is the molecular diffusivity. 
A huge amount of work has been done in the last years on 
the Kraichnan model. Due to the white-in-time charac- 
ter of the advecting velocity field, the equation for passive 
correlators of any order n are linear and closed. This al- 
lows explicit, perturbative calculations of anomalous ex- 
ponents in terms of zero-mode solutions of the closed 
equation satisfied by n-points correlation function, by 
means of developments in £ || or in 1/d Q , with d the 
physical space dimensionality. 

The connection between anomalous scaling and zero 
modes, if fascinating from one side, looks very difficult 
to be useful for the most important problem of Navier- 
Stokes eqs. In that case, being the problem non linear, 
the hierarchy of equations of motion for velocity correla- 
tors is not closed and the zero-mode approach should be 
pursued in a much less handable functional space. 
From a phenomenological point of view, a simple way 
to understand the presence of anomalous scaling is to 
think at the scalar field as made of singular scaling 



fluctuations S r 9(x) ~ r H x ) ) with a probability to de- 
velop an /i-fluctuation at scale r given by P r {h) ~ r^ h \ 
being f(h) the co-dimension of the fractal set where 
h(x) = h. This is the multifractal road to anomalous 
exponents that leads to the usual saddle-point esti- 
mate for the scaling exponents of structure functions: 
C(2n) = min?,, (2nh + f(h)) ||. In this framework, high 
order structure functions are dominated by the most in- 
tense events, i.e. fluctuation characterized by an expo- 
nent hrain'- lfm, WO o C( n ) = n hmin- The emergence of 
singular fluctuations, at the basis of the multifractal in- 
terpretation, naturally suggests that instantonic calculus 
can be used to study such special configurations in the 
system. Recently, instantons have been successfully ap- 
plied in the Kraichnan model to estimate the behaviour 
of high-order structure functions when d(2 — £) 3> 1 ||, 
and to estimate PDF tails for £ = 2 [|| . 
In this letter, we propose an application of the instan- 
tonic approach in random shell models for passive scalar 
advection, where explicit calculation of the singular co- 
herent structures can be performed. 
Let us briefly summarized our strategy and our main 
findings. First, we restrict our hunt for instantons to 
coupled, self-similar, configurations of noise and passive 
scalar, a plausible assumption in view of the multifractal 
picture described above. We develop a method for com- 
puting in a numerical but exact way such configurations 
of optimal Gaussian weight for any scaling exponent h. 
We find that h cannot go below some finite threshold 
h m in(£)- We compare h m in(£) a t varying £ given from 
the instantonic calculus with those extracted from nu- 
merical simulation showing that the agreement is perfect 
and therefore supporting the idea that self-similar struc- 
tures gouvern high-order intermittency. 
Second, assuming that these localized pulse-like instan- 
tons constitute the elementary bricks of intermittency 
also for finite-order moments we compute their dressing 
by quadratic fluctuations. We obtain in this way the first 
two terms of the function f(h) via a "semi-classical" ex- 
pansion. Let us notice that a rigorous application of the 
semi-classical analysis would demand for a small param- 
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eter controlling the rate of convergence of the expansion, 
like 1/n where n is the order of the moment [0] or 1/d, 
where d is the physical space dimension [^J. As we do 
not dispose of such small parameter in our problem, the 
reliability of our results concerning the statistical weight 
of the /i-pulses can only be checked from an a posteriori 
comparison with numerical data existing in literature. 
At the end of this communication, we will present some 
preliminary results on such important issue, while much 
more extensive work will be reported elsewhere. 
Shell models are simplified dynamical models which have 
demonstrated in the past to be able to reproduce many 
of the most important features of both velocity and pas- 
sive turbulent cascades ||. 

The model we are going to use is defined as follows. 
First, a shell-discretization of the Fourier space in a 
set of wavenumbers defined on a geometric progression 
k m = ko\ m is introduced. Then, passive increments at 
scale 7*771 ~~ kjy. 1 are described by a real variable m (t). 
The time evolution is obtained according to the follow- 
ing criteria : (i) the linear term is purely diffusive and 
is given by —nk m 9 m ; (ii) the advection term is a com- 
bination of the form k m 9 m >u m » , where u m are random 
Gaussian and white-in-time shell-velocity fields; (iii) in- 
teracting shells are restricted to nearest-neighbors of m; 
(iv) in the absence of forcing and damping, the model 
conserves the volume in the phase-space and the energy 
E = J2 m \0m\ 2 - Properties (i), (ii) and and (iv) are valid 
also for the original equation (pi) in the Fourier space, 
while property (iii) is an assumption of locality of in- 
teractions among modes, which is rather well founded 
as long as <C £ <C 2. The simplest model exhibiting 
inertial-range intermittency is defined by : 

[-T + nk m ] ^ (t) = c m 6 m -i(t)u m (t) + 



(t)u m _i(i) + 6 lm (fr(t), 



(2) 



with 0,77 = — Cm— i = k m , and where the forcing term 
acts only on the first shell. Following Kraichnan, we also 
assume that the forcing term <f>{t) and the velocity vari- 
ables u m (t) are independent Gaussian and white-in-time 
random variables, with the following scaling prescription 
for the advecting field: 



(u m (t)u n {t')) = 5(t-t')k- £ S r , 



(3) 



Shell models have been proved analytically and non- 
perturbatively |J to possess anomalous zero modes sim- 
ilarly to the original Kraichnan model ([!]). 
The role played by fluctuations with local exponent h(x) 
in the original physical space model is here replaced by 
the formation at larger scale of structures propagating 
self-similarly towards smaller scales. The existence in 
the inviscid unforced problem of such solutions associ- 
ated with the appearance of finite time singularities is a 
The analytical resolution of the instantonic problem even 
in the simplified case of shell models is a hard task. In 
iDL a numerical method to select self-similar instantons 



in the case of a shell model for turbulence, has been intro- 
duced. In the following, we are going to apply a similar 
method to our case. 

We rewrite model (0) in a more concise form: 



d6 
~dt 



M[b]0. 



(4) 



The scalar and velocity gradient vectors, and b, are 
made from the variables 6*777 and k m u m . As far as inertial 
scaling is concerned, we expect that some strong uni- 
versality properties apply with respect to the large scale 
forcing. Indeed, forcing changes only the probability with 
which a pulse appears at large scale, but not its inertial 
range scaling behaviour, Pk m (h) ~ km* • So, as we are 
interested only in the evaluation of f(h), we drop the 
forcing and dissipation in (|j). The matrix M[b] is linear 
in b and can be obviously deduced from (|J) . The stochas- 
tic multiplicative equation (|j) must be interpreted a la 
Stratonovich. Nevertheless, once the Ito-prescription for 
time discretization is adopted, the dynamics gets Marko- 
vian and a path integral formulation can then be easily 
implemented. This changes (0) into: 



dO 
~~dl 



5D0 + M[b]0, 



(5) 



where D is a diagonal matrix (Ito-drift) D mr , 
and B is a positive constant. 

As we said before, we are looking for coherent structures 
developing a scaling law 6 m ~ k m h as they propagate to- 
wards small scales in the presence of a velocity realization 
of optimal Gaussian weight. The probability to go from 
one point to another in configuration space (spanned by 
0) between times U and tf can be written quite gener- 
ally as a path integral over the three fields b, 0, p of 

_<?rhfl 1 ~~ f tf C[h,9,p]dt 

the exponential e a \. D '°>Pi — e Jh ^ where the 

Lagrangian C is given by the equation: 



1 dO 
£(b, 9,p) = -b.D-^ + p.(— + BB9 - M[b]0), (6) 

and p is an auxiliary field conjugated to which enforces 
the equation of motion (^|). The minimization of the ef- 
fective action S leads to the following coupled equations: 



d9 

— = -5D0 + M[b]0, 
dt 

^=BDp- *M[b]p, 

with the self-consistency condition for b: 
b = D'N[0]p, 



(7) 
(8) 

(9) 



where the matrix N[0] is defined implicitly through the 
relation N[0]b = M[b]0. 

We are now able to predict the scaling dependence of 
variables b m . For a truly self-similar propagation, the 
cost in action per each step along the cascade must be 
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constant. The characteristic turn-over time required by 
a pulse localized on the m — th shell to move to the next 
one can be dimensionally estimated as l/(u m k m ) = b m l . 
Recalling the scaling dependence of D and the definition 



Cdt 



-(2-0, 



of action (|6|), we expect: AS — 

We can thus deduce that b m ~ k^^. 
Let us now discuss how to explicitly find solutions of 
the above system of equations. Clearly, there is no hope 
to analytically find the exact solutions of these deter- 
ministic non linear coupled equations. Also numerically, 
the problem is quite delicate, because (0) and (||) are 
obviously dual of each other and have opposite dynami- 
cal stability properties. This phenomenon can be hardly 
captured by a direct time integration. To overcome this 
obstacle, in (lj it has been proposed a general alterna- 
tive scheme which adopts an iterative procedure. For a 
given configuration of the noise, each step consists in inte- 
grating the dynamics of the passive scalar (0) forward in 
time to let emerge the solution of optimal growth. Con- 
versely, the dual dynamics of the auxiliary field (||) is 
integrated backward in time, along the direction of min- 
imal growth in agreement with the prediction deduced 
from (||): ||p|| ~ Then the noise b can be re- 

computed by the self-consistency equation (^) and the 
process is repeated until the convergence is reached. 
Self-similar passive solutions must be triggered by self- 
similar noise configuration: 



b m (t) 



1 



(t* - t) 



F(k m ^(t*-t)), 



(10) 



where t* is the critical time at which a self-similar so- 
lution reaches infinitesimally small scales in absence of 
dissipation. To overcome the non-homogeneity of time 
evolution seen by these accelerating pulses, we introduce 
a new time variable r = — logft* — t) . Then, the advert- 
ing self similar velocity field (JTG) can be rewritten under 
the form: b(r) = e T C(r) where C m (r) is still the veloc- 
ity gradient field, but expressed in a different time scale, 
such that: C m (r) = F(m (2 - £) log A - r). 
The sought self-similar solutions appear in this represen- 
tation as traveling waves, whose period T = (2 — £) log A 
is fixed by the scaling consideration reported above. In 
this way, we can limit the search of solutions on the time 
interval [0 — T] , and the action at the final time tf = raT 
is deduced by S(t f ) = mS(T). 

Then comes the main point of our algorithm. For a fixed 
noise configuration C, the field 6 must be the eigenvector 
associated to the maximal (in absolute value) Lyapunov 
exponent a max of the Floquet evolution operator: 

U(T;0) = T_i exp f {-BD e - T + M[C(r)])dr. (ff) 
Jo 

Here TL% denotes the translation operator by one unit 
to the left along the lattice. Similarly, the auxiliary field 
must be the eigenvector associated with the Lyapunov 
exponent —a max of the inverse dual operator t U^ 1 . 



Starting from an initial arbitrary traveling wave shape 
for C(t) with period T, we have computed the passive 
scalar and its conjugate fields at any time between and 
T, by diagonalization of operator U, recomputed the ve- 
locity gradient field C(t) from the self-consistency equa- 
tion (||) and iterated this procedure until an asymptotic 
stable state, 9°, p°, C°, was reached. The scaling expo- 
nent 9 m ~ k m h for the passive scalar can be deduced by 
e a m {h) ~ e m,T — T , so that h=(£- 2)a max . Note that h 
is bound to be positive due to the conservation of energy. 
In our algorithm, the norm of the gradient velocity field 
C(0) acts as the unique control parameter in a one to one 
correspondence with h. The action S '(h) is, in multifrac- 
tal language, nothing but the first estimate of f(h) curve 
based only on the contribution of all pulse-like solutions, 
more precisely f(h) = S°(h)/ln\. 

We now turn to the presentation and discussion of our 
main result. By varying the control parameter, we ob- 
tain a continuum of exponents in the range h m i n (^) < 
h < h max (£,). The simple analysis of the /i-spectrum 
allows predictions only for observable which do not de- 
pend on the f(h) curve, i.e. only on the scaling of T 2n 
for n — > oo, (C(n) ~ h m i n n for n large enough). 
Unfortunately, high order exponents are the most dif- 
ficult quantities to be extracted from numerical or ex- 
perimental data. Nevertheless, thanks to the extreme 
simplicity of shell models, very accurate numerical simu- 
lations have been done |L2| at different values of £ and in 
some cases a safe upper bound prediction on the asymp- 
totic of exponents could be extracted. 




FIG. 1. Behavior of the asymptotic slope h m i n as function 
of 7 = 2 — £ : (+) results from the instantonic calculus for 
the model of Ref. Q , (*) upper bound value for h min as 
extracted from numerical integration of the same model. 

To compare our results with the numerical data existing 
in literature, we have analyzed the shell-model version of 
passive advection proposed in |12|. In Fig.f , we show the 
hmin curve obtained at various £ from instantonic cal- 
culation, together with the h™™ values extracted from 
direct numerical simulation of the quoted model Jl^] per- 
formed at two different values of £: the agreement is 
good. Our calculation predicts, within numerical errors, 
the existence of a critical £ c ~ 1.75 above which the min- 
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imal exponent reaches the lowest bound h m i n = 0. 
This goes under the name of saturation and it is the sig- 
nature of the presence of discontinuous-like solutions in 
the physical space S r 6 ~ r°. Theoretical || and numeri- 
cal |l3| results suggest the existence of such effect in the 
Kraichnan model for any value of £. The existence of sat- 
uration in this last is due to typical real-space effects and 
therefore it is not surprising that there in not a complete 
quantitative analogy with the shell-model case. 
Let us now present the other -preliminary- result, i.e. the 
role played by instantons for finite-order structure func- 
tions. If we just keep the zero-th order approximation for 
f(h) = So (h)/ log A, we get the ( n curve shown in Fig. 2, 
which is quite far from the numerical results of (the 
asymptotic linear behavior is in fact not even reached in 
the range of n represented on the figure). In order to 
get a better assessment of the true statistical weight of 
the optimal solutions, we computed the next to leading 
order term in a "semi-classical" expansion. Fluctuations 
around the action were developed to quadratic order with 
respect to 9°, p°, C°, and the summation over all per- 
turbed trajectories leading to the same effective scaling 
exponent for the 6 field after m cascade steps was per- 
formed. It turns out (see Fig. 2) that the contribution in 
the action of quadratic fluctuations, Si(h), greatly im- 
proves the evaluation of £( n )- 



also on not too intense fluctuations. 
In conclusion, we have presented an application of the 
semi-classical approach in the framework of shell models 
for random advection of passive scalar fields. Instantons 
are calculated through a numerically assisted method 
solving the equations coming from probability extrema: 
the algorithm has revealed capable to pick up those con- 
figurations giving the main contributions to high order 
moments. Of course, we are far from having a systematic, 
under analytical control, approach to calculate anoma- 
lous exponents in this class of models. Nevertheless, the 
encouraging results here presented raise some relevant 
questions which go well beyond the realm of shell-models. 
To quote just one, we still lack a full comprehension of 
the connection between the usual multiplicative-random 
process and the instantonic approaches to multifractality: 
in particular, it is not clear what would be the prediction 
for multi-scale and multi-time correlations of the kind 
discussed in fll5f within the instantonic formulation. 
It is a pleasure to thank J-L. Gilson and P. Muratore- 
Ginanneschi for many useful discussions on the sub- 
ject. LB has been partially supported by INFM (PRA- 
TURBO) and by the EU contract FMRX-CT98-0175. 
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FIG. 2. Exponents ( n of structure functions of order n vs 
n for £ = 1: (*) data from direct numerical simulation of 
Ref. [M; (+) Legendre transform of the action f(h) — So', 
(•) Legendre transform of the action f(h) = So + Si. The 
solid line corresponds to the dimensional scaling. 

Naturally, in the absence of any small parameter in the 
problem, we cannot take for granted that the next correc- 
tion^) would not spoil this rather nice agreement with 
numerical data. But the surprising fact that Sq + Si is 
strongly reduced with respect to So, even for the most 
intense events, does not imply by itself a lack of consis- 
tency of our computation. In any case, the prediction of 
the asymptotic slope of the curve, based on the value 
hmin is obviously valid beyond all orders of perturbation. 
Moreover, for values of £ > 1, we find that the second or- 
der exponent extracted from our calculation is in good 
agreement the exact result Q 2 — 2 — £, suggesting that our 
approach is able to give relevant statistical information 
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